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We analyze the electronic Green's functions in the superconducting ground state of the t—J model 
using Gutzwiller-projected wave functions, and compare them to the conventional BCS form. Some 
of the properties of the BCS state are preserved by the projection: the total spectral weight is 
continuous around the quasiparticle node and approximately constant along the Fermi surface. On 
the other hand, the overall spectral weight is reduced by the projection with a momentum-dependent 
renormalization, and the projection produces electron-hole asymmetry in renormalization of the 
electron and hole spectral weights. The latter asymmetry leads to the bending of the effective Fermi 
surface which we define as the locus of equal electron and hole spectral weight. 
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I. INTRODUCTION 

Shortly after the discovery of superconductivity 
in copper-oxyde compounds^ Anderson proposed a 
Gutzwiller-projected BCS wave function which would 
describe the superconducting ground state of high- 
temperature superconductors^ The variational approach 
to superconducting cuprates based on Anderson's origi- 
nal proposal has since had a lot of success, while the 
strong Coulomb repulsion and the non-perturbative na- 
ture of the problem make other approaches extremely dif- 
ficult. Interest in projected wave functions as variational 
ground states for cuprate superconductors was initiated 
by several research groups in late 80's^ and lead to a 
considerable activity in the field. The projected wave 
functions show large overlap with exact ground states on 
small clusters and have low variational energies for the 
t—J model£& Furthermore, several experimental prop- 
erties of cuprates like the zero-temperature phase dia- 
gram and d-wave pairing symmetry are extremely well 
predicted within this approach^ 

Due to considerable progress of the experimental tech- 
nique of Angle Resolved Photoemission Spectroscopy 
(ARPES) on cuprates, more and more high-quality 
data on the low-lying spectral properties of these com- 
pounds have been made available in recent years^ Ex- 
perimentally, the low-energy excitations of superconduct- 
ing cuprates are known to resemble BCS quasiparticles 
(QPs). 10 It is therefore interesting to theoretically ex- 
plore the wave function of projected QP excitations and 
compare them to unprojected BCS QPs. The most ap- 
parent differences are the the doping dependency of the 
nodal Fermi velocity and the renormalizations of the 
nodal QP spectral weight and of the current carried by 
Qp S( 7 i ii i i2 4 i3 i i4 j n .j. ne p resen ^ paper we further analyze 

the properties of the superconducting ground state and 
the QP excitations with the Variational Monte Carlo 
technique (VMC)pi^ We calculate the equal-time Green's 
functions, both normal and anomalous, in the Gutzwiller- 
projected state and derive from them the QP spectral 
weights for addition and removal of an electron at zero 
temperature. The main conclusion of our study is that, 



due to a non-trivial interplay of superconductivity and 
strong Coulomb repulsion (projection), the electron and 
hole spectral weights are renormalized differently. A nat- 
ural way to describe this asymmetry is to define the "ef- 
fective Fermi surface" as the locus of points where the 
electron and the hole spectral weights are equal. Thus de- 
fined Fermi surface acquires an additional outward bend- 
ing in the anti-nodal region as compared to the original 
unprojected Fermi surface. This bending is a signature of 
a deviation from the BCS theory and may be responsible 
for the shape of the Fermi surface observed in ARPES ex- 
periments. The validity of Luttinger's rulei& in strongly 
interacting and superconducting materials has been ques- 
tioned experimentally and theoretically recently. 917 - 18 
Our findings provide further indication of its inapplica- 
bility in strongly correlated superconductors. 

The paper is organized as follows. Section ITT1 contains 
the definition of the model and wave functions used in our 
calculations. In Section IlIII we first present some exact 
relations for projected wave functions, we then describe 
our results on the QP spectral weights. Section \N\ is 
devoted to the calculation of the equal-time anomalous 
Green's function. Finally, Section E] defines the "effec- 
tive Fermi surface" and discusses its deviation from the 
unprojected one. 



II. THE MODEL 

In the tight-binding description, the cuprates are mod- 
eled by electrons hopping on a square lattice. The ap- 
propriate model is the t—J Hamiltonian: 

H t .j = -t £ P G clc ja P G + jJ2(Si-S j -^) 

<i,j>,o <i,j> 

(1) 

acting in the Hilbert-space with less than two electrons 
per site. Here n ia = c\ a c il7 , Si = ^^cTao'Cia' , c\ a is the 
electron creation operator in the Wannier state at site i, 
and er are the Pauli matrices. The no-double occupancy 
is preserved by the Gutzwiller projector Pq = FL;(1 — 
n lT ni|). 
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The t—J model can be viewed as the large U limit of the 
one-band Hubbard model, neglecting the 3-site-hopping 
term. Provided that the model is analytic in t/U, dou- 
bly occupied sites can be re-introduced perturbatively to 
recover the full Hilbert space of the Hubbard modelAiS 
Although the inclusion of these corrections present no 
major difficulty, we chose to neglect them here. In most 
quantities, only small corrections arise from finite double 
occupancyr*^ which makes this approach to the large- 
U Hubbard model consistent. Furthermore, it has been 
argued that the t-J model is in fact more appropriate 
than the one-band Hubbard model in describing the CuO 
planes 

We consider the usual variational ground stated 



\V H ) = P H PG\dBCS(A, f j,)) 



(2) 



where Pr is the particle number projector on the sub- 
space of L 2 — H electrons; L 2 is the total number of 
sites. Both hole number and number of sites are even. 
\dBCS) is the ground state of the BCS mean field Hamil- 
tonian with nearest neighbor hopping and d-wave pair- 
ing symmetry on the square lattice: \dBCS) = Tl k (u k + 
v fc c feT c -fcl)l°) K n fc-ff7fc<r|0). 7fc<r = u k c klJ + <Tv k c f _ ksl 

< = l(i - 4 = l(i + &)> E h = v/^Tai, a = 

— 2(cos(fe c ) + coa(ky)) — fj,, A k = A(cos(k x ) — cos(k y )). 
The wave function (0 has two free parameters: A and /i. 
These variational parameters are chosen to minimize the 
energy of the t—J Hamiltonian Q for the experimentally 
relevant value J/t = 0.3 and for every doping level. We 
use the optimized parameters and the cluster geometry 
from Ref. M 

The following ansatz is used for the excited 
statesiWV 2 ! 13 ! 14 



\^H,k,a) = P H P G lL\dBCS) 



(3) 



In the following, the normalized versions of J2J and © 
will be denoted by \H) and \H,k,a), respectively: 

\H) = H^ffir 1 !^) (4a) 
\H,k,a) = I^H^ir 1 !*^) (4b) 

The many-particle wave function (J2J) , sometimes called 
Anderson's Resonant Valence Bond (RVB) wave func- 
tion, implements both strong electron correlations and 
superconductivity. It is known to have a considerable 
overlap with the true ground state of the t-J model at 
non-zero hole doping on small clusters^SiSLS There is 
also numerical support from exact diagonalization studies 
indicating well defined BCS-like QPs as low-energy exci- 
tations of the t-J models Therefore, the excited trial 
states © are expected to be close to the true excitations 
of the t—J model. However, here we are more interested 
in the physical content of the proposed wave functions 
than in their closeness to the eigenstates of a particular 
Hamiltonian. 37 



III. QUASIPARTICLE SPECTRAL WEIGHTS 

The QP spectral weights are defined as the overlap 
between the bare electron/hole added ground state and 
the QP excitations of the model: 



Zi H = \(HTl,k,a\c±\H) 



(5) 



where c fc are the bare electron creation/anihilation op- 
erators. 

It is well-known that the QP spectral weight for adding 
an electron (equal-time normal Green's function) can be 
calculated from the ground state spectral functioniii 



7+ - 



1 + X 



( c L c ka) 



(6) 



where x = H/L 2 is the hole concentration. 

The QP spectral weight for adding a hole is more dif- 
ficult to calculate. It is useful to note that it can also be 
calculated from ground state expectation values, 



Z k,H+l Z kM-l = I $ mH 



(7) 



where ^>k,H is the superconducting order parameter 
(equal-time anomalous Green's function) 



®k,H = (H + \\c k]C - ki \H -I) 



(8) 



Relation (JJJ can be proven by algebraic manipulations 
with the Gutzwiller projector. It is exact in finite systems 
and is also valid in the thermodynamic limit: 



ztz k = |$ fc | 



(9) 



Remarkably, this relation holds for both the unprojected 
BCS state (with Z+ = \v k \ 2 , Z k = \u k \ 2 ) and the fully 
projected wave function. It would be interesting to ex- 
plore if this relation is valid in a more general case of 
superconducting systems or if it is just a peculiarity of a 
certain class of wave functions. 

Further, we define the total spectral weight 



rftot 



Z k + z k 



(10) 



The main contribution to Z k ot is given by Z k outside the 
Fermi surface and by Z k inside. We can prove an exact 
upper bound on Z k ot : 



< 



1 + x 



(11) 



A proof may be performed by defining the two states 

\a) = c feT |ff-l) (12a) 

|6> = P G cl hL \H+l) (12b) 

(with this definition, \b) is proportional to \H, k, f)). Us- 
ing |[BJ and © , we show that 



/(o|o> (a\b) 
\(b\a) (b\b) 



(13) 



^ - zt] zi |$ fc | 2 = z£ f 1 -^ zr 
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On the other hand, the same determinant equals 



l\M®\b) 



\b)®\a)\\ 2 >0, 



(14) 



which completes the proof. 

Numerically, we compute the spectral weight H1 
by first computing Z\ 



k,H+l 



and <&k,Hi and then using 
(J7J). The disadvantage of this method is large error bars 
around the center of the Brillouin zone where both 
and $ fc are small (recently performed calculations of Z k 
by direct sampling of the excited states are free from 
this problem)^ However, our precision is sufficient to 
establish that the total spectral weight is a smooth 
function and has no singularity at the nodal point. 

Technically, &k,H is computed as \&k,H 2 = &t H^k h > 
where 



3> 



k.H 



3> 



H+l\ck]C- k j\^ H-l) 

ll*H+l|| 2 
(*H + l|Cfe-fC_fe4_|*H-_i) 



k.H 



H-l 



(15a) 
(15b) 



Both matrix elements can be computed within the usual 
Metropolis algorithm^ 



In Figs . ^ and [21 we plot the spectral weights Z~£ , Z^ , 
and Z^* along the contour — > (0, tt) — » (tt, tt) — > 
in the Brillouin zone for different doping levels. Fig- 
ure shows the contour plots of Zt' 1 in the region of the 
Brillouin where our method produces small error bars^ 
From these data, we can make the following observations: 

• In the case of an unprojected BCS wave function, 
the total spectral weight is constant and unity over 
the Brillouin zone. Introducing the projection op- 
erator, we see that for low doping (x ~ 3%), the 
spectral weight is reduced by a factor up to 20. The 
renormalization is asymmetric in the sense that the 
electronic spectral weight Z£ is more reduced than 
the hole spectral weight Z^ . For higher doping 
(x ~ 23%), the spectral weight reduction is much 
smaller and the electron-hole asymmetry decreases. 

• Since there is no electron-hole mixing along the 
zone diagonal, the spectral weights Z^ and Z^ 
have a discontinuity at the nodal point. Our data 
shows that the total spectral weight is continuous 
across the nodal point. Strong correlations does 
not affect this feature of uncorrelated BCS the- 
ory. Recently, it has been argued in Ref. |2(| that 
the total spectral weight of the projected (non- 
superconducting) Fermi sea should be continuous 
across the Fermi surface. This is consistent with 
our result. 

The intensities measured in ARPES experiments are 
proportional to the spectral weights of the low-energy 
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FIG. 1: (Color online) QP spectral weights for 6 holes (upper 
plot, x ~ 3%) and 22 holes (lower plot, x ~ 11%) on 196 
sites. The spectral weights are plotted along the contour — > 
(0, tt) — > (7T, tt) — > (shown in inset). Plus signs (+, blue 
online) denote the spectral weight Z£ , crosses (x) denote 
Z^ , error bars are shown. Solid dots (black) denote their 
sum, the total spectral weight Z^, error bars not shown. On 
the horizontal axis, the star (*) denotes the intersection with 
the unprojected Fermi surface along the — > (0, tt) direction; 
the thick dot is the nodal point. Z^ and Z^ jump at the 
nodal point, while Z£ ot is continuous. The intersection with 
the effective Fermi surface (see section [VJ is marked by an 
arrowhead. On the diagonal (last segment), k is given in 
units of y/2. 



0.5 
0.4 
0.3 
0.2 
0.1 



1 1 — 




T 1 


1 1 1 1 


1 1 1 


1 










-t- - 








« +" 
-+- 


+ + + ^ + + + + 














-+- 










1 


a & a a a * a mj^J^itJ^± 



.2 .4 



1 .2 .4 .6 .8 1 .8 .6 .4 .2 



0.4 
0.3 
0.2 
0.1 



1 — ■ 


— 1 


-\ 1 

4 " * ' • 

: -+- 


1 1 1 1 


— i — i — i — 

++++ 


i 














■++^+ 


-f- 


. * *«# * * m * * 




n— n- 



.2 .4 .6 



1 .2 .4 .6 .8 1 

k/n 



.6 .4 .2 



FIG. 2: (Color online) Same plot as Fig.^of the QP spectral 
weights for 34 holes (upper plot, x ~ 17%) and 46 holes (lower 
plot, x ~ 23%) on 196 sites. 
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QPs£ In Ref. IToL the spectral weights of a slightly over- 
doped sample of Bi2223 were measured along the cut 
(-7T, 0) — > (tt/2, 7r/2); an almost constant total spectral 
weight was reported in this experiment. It can be seen 
from Fig. [5] that the total spectral weight is approxi- 
mately constant along this cut, so the experimental result 
agrees with this property of projected wave functions. 

An anisotropy of the ARPES intensity along the ex- 
perimental FS (the so-called nodal-antinodal dichotomy) 
has been reported in a series of experiments^^ Experi- 
mentally, the spectral weight measured in the anti-nodal 
region is suppressed in undcrdoped compounds, while it 
is large in the optimally doped and overdoped region. 
Usually, this effect is associated with formation of some 
charge or spin order, static or fluctuating one. From 
Fig. [5] we see that a similar (but much weaker) tendency 
can be observed in the framework of Gutzwillcr-projected 
wave functions. The experimentally observed effect is 
much stronger and a claim that the nodal-antinodal di- 
chotomy can be explained within this framework would 
be too hasty. 



IV. SUPERCONDUCTING ORDER 
PARAMETER 
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FIG. 3: (Color online) Doping dependency of the nearest- 
neighbor superconducting order parameter $jj (calculated in 
the 14x14 system). The error bars are smaller than the sym- 
bol size. The same quantity calculated in the Gutzwiller ap- 
proximation is also shown for comparison. The variational 
parameter A is shown with the scale on the right. 

In Fig. 13 we plot the nearest-neighbor superconduct- 
ing correlation <E>y (the Fourier transform of defined 
in (JSJ) as a function of doping. This curve shows close 
quantitative agreement with the result of Ref. 0, where 
the authors extracted the superconducting order parame- 
ter from the long range asymptotics of the nearest neigh- 



bor pairing correlator, lim r ^ 00 (coC5cJc^ +(5 ). With the 
method employed here, we find the same qualitative and 
quantitative conclusions of previous authorsi&i vanish- 
ing of superconductivity at half filling x — > and at the 
superconducting transition on the overdoped side x c ~ 
0.3i2& The optimal doping is near x op t — 0.18. In the 
same plot we also show the commonly used Gutzwiller 
approximation where the BCS order parameter is renor- 
malized by the factor g t = £ The Gutzwiller approxi- 
mation underestimates the exact value by approximately 
25%. 
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FIG. 4: (Color online) Pairing correlation <E>fc in the Brillouin 
Zone at differen doping levels: x ~ 3% (upper left), 11% (up- 
per right), 17% (lower left) and 23% (lower right) in a 14 x 14 
system with periodic-antiperiodic boundary conditions. The 
dashed line is the unprojected FS. 

In Fig.0] we show contour plots of the superconducting 
order parameter for four values of doping. It resem- 
bles qualitatively the unprojected d-wave pairing ampli- 
tude, but is somewhat distorted due to the particle-hole 
asymmetry (see discussion in the previous and the fol- 
lowing sections). 



V. FERMI SURFACE 

In strongly interacting Fermi systems, the notion of a 
Fermi surface (FS) is not at all clear. There are however 
several experimental definitions of the FS. Most com- 
monly, kp is determined in ARPES experiments as the 
maximum of |Vfcnfe| or the locus of minimal gap along 
some cut in the As-planed The theoretically better defined 
locus of rife = 1/2 is also sometimes used. The various 
definitions of the FS usually agree within the experimen- 
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tal uncertainties. Recently, the different definitions of the 
FS were theoretically analyzed in Refs. Il7l and fTst 

In our present work, we propose an alternative defi- 
nition of the Fermi surface based on the ground state 
equal-time Green's functions. In the unprojected BCS 
state, the underlying FS is determined by the condition 
\uk\ 2 — \vk\ 2 - We will refer to this as the unprojected FS. 
Since \uk\ 2 and \vk\ 2 are the residues of the QP poles 
in the BCS theory, it is natural to replace them in the 
interacting case by Z£ and , respectively. We will 
therefore define the effective FS as the locus Z£ = Z^ . 




k In kin 



FIG. 5: (Color online) Contour plots of the total QP spectral 
weight Zl ot . The effective FS (dash line) and unprojected 
FS (dash-dot line) are also shown. The doping levels are 
x ~ 3% (top left), 11% (top right), 17% (bottom left) and 
23% (bottom right). The + signs indicate points where the 
values are known within small error bars. 

In Fig. [3]we plot the unprojected and the effective FS 
which we obtained from VMC calculations. The contour 
plot of the total QP weight is also shown. It is interesting 
to note the following points: 

• In the underdoped region, the effective FS is open 
and bent outwards (hole- like FS). In the overdoped 
region, the effective FS closes and embraces more 
and more the unprojected one with increasing dop- 
ing (electron- like FS). 

• Luttinger's rule^ is clearly violated in the under- 
doped region, i.e. the area enclosed by the effective 
FS is not conserved by the interaction; it is larger 
than that of the unprojected one. 

• In the optimally doped and overdoped region, the 
total spectral weight is approximately constant 



along the effective FS within error bars. In the 
highly underdoped region, we observe a small con- 
centration of the spectral weight around the nodal 
point (~ 20%). 

A large "hole-like" FS in underdoped cuprates has 
also been reported in ARPES experiments by several 
groupsi 27 i 29 i 30 

It should be noted that a negative next-nearest hop- 
ping t' would lead to a similar FS curvature as we find in 
the underdoped region. We would like to emphasize that 
our original t-J Hamiltonian as well as the variational 
states do not contain any t' . Our results show that the 
outward curvature of the FS is due to strong Coulomb 
repulsion, without need of t'. The next-nearest hopping 
terms in the microscopic description of the cuprates may 
not be necessary to explain the FS topology found in 
ARPES experiments. Remarkably, if the next-nearest 
hopping t' is included in the variational ansatz (and not 
in the original t-J Hamiltonian), a finite and negative 
t' is generated, as it was shown in Ref. 0. Apparently, 
in this case the unprojected FS has the tendency to ad- 
just to the effective FS. A similar bending of the FS was 
also reported in the recent analysis of the current car- 
ried by Gutzwiller- projected QPs. 14 A high-temperature 
expansion of the momentum distribution function nj- of 
the t-J model was done in Ref . |32j where the authors find 
a violation of Luttinger's rule and a negative curvature 
of the FS. Our findings provide further evidence in this 
direction. 

A natural question is the role of superconductivity 
in the unconventional bending of the FS. In the limit 
A — > 0, the variational states are Gutzwiller-projected 
excitations of the Fermi sea and the spectral weights 
are step-functions at the (unprojected) FS. In a recent 
paper— it was shown that hm.. .+ Z^ = lim fe - Z^ 
for the projected Fermi-sea, which means that the un- 
projected and the effective FS coincide in that case. This 
suggests that the "hole-like" FS results from a non-trivial 
interplay between strong correlation and superconduc- 
tivity. At the moment, we lack a qualitative explana- 
tion of this effect, however it may be a consequence of 
the proximity of the system to the non-superconducting 
"staggercd-flux" stat o 21 ! 33 or to anti-ferromagnetismLM 
near half-filling. 
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